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We investigate the viability of the Galactic halo model of 7-ray 
. bursts by calculating the spatial distribution of neutron stars born 

with high velocities in the Galactic disk, and comparing the result- 
ing brightness and angular distribution with the BATSE data. We 
find that the Galactic halo model can reproduce the BATSE peak 
bJQ' flux and angular distribution data for neutron star kick velocities 

^ , ^ 800 km s^^, source turn-on ages ^ 10 Myrs, and sampling depths 

<^ ■ 100 kpc ^ dmax ^ 400 kpc. 
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INTRODUCTION 



o 



Gamma-ray bursts (GRBs) continue to confound astrophysicists nearly a 
quarter century after their discovery (^. Before the launch of CGRO, most 
I scientists thought that GRBs came from magnetic neutron stars residing in a 
^ thick disk (having a scale height of up to ~ 2 kpc) in the Milky Way (H,^). The 
I data gathered by BATSE showed the existence of a rollover in the cumulative 
I brightness distribution of GRBs and that the sky distribution of even faint 
. GRBs is consistent with isotropy (^,||). This rules out a thick Galactic disk 
Qh' source population. 

Q . Consequently, the primary impact of the BATSE results has been to inten- 
■ sify debate about whether the bursts are Galactic or cosmological in origin. 
c/3 Galactic models attribute the bursts primarily to high-velocity neutron stars 
in an extended Galactic halo, which must reach one fourth or more of the dis- 
J> tance to M31 {dusi ~ 690 kpc) in order to avoid any discernible anisotropy 
(^,0). Cosmological models place the GRB sources at distances d ^ 1 — 3 
Gpc, corresponding to redshifts z ~ 0.3 — 1. A source population at such 
d [ large distances naturally produces an isotropic distribution of bursts on the 
sky, and the expansion of the universe or source evolution can reproduce the 
observed rollover in the cumulative brightness distribution (^. 

Recent studies (HJl^) have revolutionized our understanding of the birth 
velocities of radio pulsars. They show that a substantial fraction of neutron 
stars have velocities that are high enough to produce an extended halo around 
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the Milky Way like that required by Galactic halo models of GRBs (11). 

Podsiadlowski, Rees, and Ruderman ( p^ have carried out pioneering cal- 
culations of the spatial distribution expected for high- velocity neutron stars 
born in the Galactic disk. They consider the effects of a non-spherical halo 
potential and of M31, but neglect the effects of the Galactic disk, which we 
find is also important. 



MODELS 

We have calculated detailed models of the spatial distribution expected for a 
population of high- velocity neutron stars born in the Galactic disk and moving 
in a Galactic potential that includes the bulge, disk, and a dark matter halo. 

We use the mass distribution and potential given by Kuijken and Gilmore 
( |l^ which includes a disk, a bulge, and a dark matter halo. The densities of 
the disk and of the halo are 

-1 

(1) 

The circular velocity Vc and the Galactic disk lead to characteristic angular 
anisotropics as a function of burst brightness which provide a signature, and 
therefore a test, of high- velocity neutron star models. Prolate or oblate dark 
matter halos also produce other angular anisotropics as a function of burst 
brightness which may provide a signature of such models (p^. 

We assume that neutron stars are born with the circular velocity ~ 
220 km s~^ of the Galactic disk. Given that current knowledge of the dis- 
tribution of initial kick velocities is uncertain, we adopt a Green-function 
approach: we calculate the spatial distribution of neutron stars for a set of 
kick velocities (e.g., i^kick = 200, 400, 1400 km s~^). We follow the resulting 
orbits for up to 3 x 10^ years. 

In our initial calculations, we assume that the bursts are standard candles, 
i.e. h = S{L — Lq). We parameterize the burst-active phase by a turn-on 
age 6t and a duration At, and assume that the rate of bursting is constant 
throughout the burst-active phase. The high- velocity neutron star model then 
has four parameters: Wkick, St, At, and the BATSE sampling depth dmax. 



PD = Pd exp 
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COMPARISON BETWEEN MODELS AND DATA 



We compare the models with a carefully-selected data set that is self- 
consistent. We use only bursts that trigger on the 1024 ms timescale be- 
cause we require that all bursts lie above the counts threshold in one trigger 
timescale; the 1024 ms timescale yields the largest sampling depth, and there- 
fore imposes the strongest constraint on models, of the three BATSE trigger 
timescales. We adopt F^^'^'^, the peak flux in 1024 ms, as our measure of 
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FIG. 1. Comparison of a Galactic halo model in which neutron stars are born 
with a kick velocity of 1000 km s~^ and have a burst-active phase lasting At = 500 
million years with a carefully-selected sample of 285 bursts from the BATSE 2B 
catalogue. Panels (a) and (b) show the contours in the {5t, dmax)-plane along which 
the Galactic dipole and quadrupole moments of the model differ from those of the 
data by ± la (solid lines), ± 2(t (dashed line), and ± 3cr (short-dashed line) where 
a is the model variance; the thin line in panel (a) shows the contour where the 
dipole moment for the model equals that for the data. Panel (c) shows the contours 
in the {5t, dmax)-plane along which 32%, 5%, and 4 x 10"'^ of simulations of the 
cumulative distribution of 285 bursts drawn from the peak flux distribution of the 
model have KS deviations D larger than that of the data. Panel (d) compares the 
brightness distribution of the model shown in (a) - (c), taking 5t = 30 Myrs and 
dmax = 200 kpc, to the BATSE plus PVO data. 



burst brightness. We therefore include only bursts which have a F^^ and 
tgo > 1024 ms. We consider only bursts with -Fp^^^ > 0.35 photons cm^^ s""'^ 
in order to avoid threshold effects (p|Jl4|). We also exclude overwriting bursts, 
because the threshold is much higher for these bursts, and MAXBC bursts, 
because they have unknown positional errors. The 2B catalogue contains 285 
bursts satisfying the above criteria. This set of bursts has Galactic dipole and 
quadrupole moments (cos 6*) = 0.056±0.034, and (sin^fe-i) = -0.033±0.017, 
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FIG. 2. Sky distribution of 285 bursts drawn randomly from the angular distri- 
bution expected for the model illustrated in Figure 1 when St = 30 Myrs, and 
rfmax = 200 kpc. 



compared to the values {cos9) — —0.013, and (sin^ & — 1) = —0.005 expected 
for a uniform sky distribution, taking into account the BATSE sky exposure. 

As a first step in testing the viabihty of Galactic halo models, we have com- 
pared the Galactic dipole and quadrupole moments, {cosO) and (sin^ 6—1/3), 
of the angular distribution of bursts for the model with those for the above 
set of bursts, using x^. We have also compared the peak flux distribution for 
the model with that for the above set of bursts, using the KS test. 

These comparisons do not provide estimates of model parameters (i.e., they 
do not yield parameter confidence regions) , but are meant only to be a rough 
"goodness-of-fit" guide to models which should be tested using a more rigorous 
approach like the maximum likelihood method. 

As an illustrative example, we show in Figure 1 the results for a Galactic 
halo model in which neutron stars are born with a uniform single velocity 
1000 km s~^ and the burst-active phase has an initial burst rate r oct^ and 
lasting = 500 million years. Figure 2 shows the sky distribution of 285 
bursts drawn randomly from the angular distribution expected for the model 
with 5t = 30 Myrs, and rfmax = 200 kpc. 

Comparisons of this kind show that the high- velocity neutron star model 
can reproduce the peak flux and angular distributions of the bursts in the 
BATSE 2B catalogue for neutron star kick velocities ?;kick ^ 800 km s~^, burst 
turn-on ages St ^ 10 million years, and BATSE sampling depths 100 kpc 
^ c^max ^ 400 kpc. Moreover, comparisons of this kind show that there is 
a large region of parameter space in which these models can reproduce the 
angular distribution of the bursts in the preliminary BATSE 3B catalogue. 

In high- velocity neutron star models, the slope of the cumulative peak flux 
distribution for the brightest BATSE bursts and the PVO bursts reflects the 
space density of the relatively small fraction of burst sources in the solar 
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neighborhood. The nearness of the observed slope of the cumulative peak 
flux distribution of these bursts to -3/2, the value expected for a uniform 
spatial distribution of sources which emit bursts that are "standard candles," 
must be considered a coincidence in the high-velocity neutron star model. 
However, a spread in neutron star kick velocities, in neutron star ages at 
which bursting behavior begins, or in the burst luminosity function tends to 
produce a cumulative peak flux distribution with a slope of -3/2; beaming of 
bursts along the direction of motion of the source or evolution of the rate of 
bursting as a function of age also tends to produce a slope of -3/2. 

We find that there are many combinations of these factors which successfully 
reproduce the slope of the BATSE plus PVO peak flux distribution. For 
example, a model in which the burst luminosity function is a log normal 
distribution with a FWHM of a factor of ^ 10 and the burst-active phase has 
an abrupt ( "heaviside function" ) turn-on, one in which the kick velocities are 
distributed between 800 and 1200 km s^^ and the burst-active phase initially 
has a burst rate r = (t/St), or one in which the burst-active phase initially has 
a burst rate r = ■^{t/St)'^ work equally well. Figure Id compares the peak flux 
distribution for the last model and the BATSE-f-PVO peak flux distribution. 

M31 provides a strong constraint on the BATSE sampling distance dmax 
@. We have investigated the effects of M31 within the framework of the 
high-velocity neutron star model described above by including the distortion 
of the Galactic halo potential due to M31 and the burst sources emanating 
from M31. We find that for such models M31 imposes a limit on the BATSE 
sampling distance dmax ^ 400 kpc, even if the bursting activity of neutron 
stars lasts for more than 10^ years (p^. 
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